On 
On 

o 

Q 



Antiferromagnetic 4-d 0(4) Model 

Isabel Campos a , Luis A. Fernandez 6 and Alfonso Tarancon" 

February 1, 2008 



a) Departamento de Fisica Teorica, Facultad de Ciencias, 
Universidad de Zaragoza, 50009 Zaragoza, Spain 
e-mail: isabel, tarancon@sol.unizar.es 



£> \ b) Departamento de Fisica Teorica I, Facultad de Ciencias Fisicas, 

Universidad Complutense de Madrid, 28040 Madrid, Spain 
e-mail: laf@lattice.fis.ucm.es 

O 

Abstract 

We study the phase diagram of the four dimensional 0(4) model 
Q-i, with first and second (^2) neighbor couplings, specially in the 

(— 1 [ fa < region, where we find a line of transitions which seems to be 

second order. We also compute the critical exponents on this line at 
the point Pi = (F4 lattice) by Finite Size Scaling techniques up to a 
r> \ lattice size of 24, being these exponents different from the Mean Field 

' ones. 
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1 Introduction 



The action for the electroweak sector of the Standard Model has SU(2)xU(l) 
symmetry. If we consider the SU(2) part and take the limit in which gauge 
degrees of freedom are frozen, the resulting action is the 0(4) non linear 
a model, which has been extensively studied because its Spontaneous Sym- 
metry Breaking pattern is related to the one exhibited by SU(2) in four 
dimensions. 

The regularized version of the 0(4) a model on the lattice leads to an 
interacting continuum limit for d < 4 0, while for d > 4 the theory is 
described by free bosonic fields [Q, ||]. 

At the upper critical dimension, d — 4, deviations from Mean Field The- 
ory (MFT) are expected. The MFT predictions for the scaling of thermody- 
namic quantities are corrected by multiplicative logarithmic terms ||. || . 

Perturbatively the infrared fixed point of the Callan-Symanzik function 
/3(g) moves to the origin as the dimension becomes four ||, also the fixed 
point is now a double zero (in contrast with the d < 4 case) which is respon- 
sible for the occurrence of such logarithmic corrections. 

The existence of these corrections imply the triviality of the theory ||. 
Triviality seems to persist when gauge fields are included [0, BJ. 

The common feature of all these approaches to the so-called triviality 
problem 0, is that the self- interactions of the scalar field in the broken 
phase are weak, and they can be reasonably studied within the context of 
perturbation theory. It is generally believed that the perturbative and the 
strong coupling regime belong to the same universality class. However, for 
the non-perturbative strong coupling regime a rigorous proof is still lacking 
and we have to rely on numerical simulations ||. [TT)| 



The existence of a strongly interacting Higgs sector, with a complicated 
dynamics, rendering useless perturbation theory, is a possibility not to be 
discarded a priori. A large amount of work have been done actually in order 
to know whether or not non-perturbative effects could change the physics of 
the electro- weak symmetry breaking sector (for a review see |11|). In this 
sense, concerning the universality class of the RP^" 1 models, the role of 



non-perturbative effects needs to be clarified ||T3 . 

Antiferromagnetism (AF) has been considered in a great variety of models 
in order to find properties not present in the purely ferromagnetic (FM) 
systems fL2 , |14[] In the context of High T c Superconductivity, AF seems to 
play an essential role. The transition from paramagnetic to non purely FM 
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I§ 0, 0,0. 



ordered phases has been studied in two dimensional models 

In four dimensions, in diluted systems recently new critical exponents 
have been obtained 
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Also in d = 4 a previous study of the AF Ising 
shows the existence of an AF phase non trivially equivalent to the 
standard FM one. However no new critical behavior was evidenced in this 
work. 

Also in four dimensions, competing interactions have been considered in 
order to study the multicritical point of the Yukawa models. At this multicrit- 
ical point four phases meet (FM, AF, Ferrimagnetic and PM). The question 
of whether or not it would be possible to define a non-trivial continuum limit 
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at this point remains still an open problem |2T|, | 

It is not clear the role that AF can play in the formulation of QFT, 
nevertheless it is worthwhile a careful study of this kind of models since they 
are known to have very rich phase diagrams, and presumably new universality 
classes could appear in which alternative formulations of continuum QFT 
should be possible. However, when defining a theory with AF couplings 
one has to be aware of the fact that higher order derivatives tends to violate 



reflection positivity [g5j, . A possibility is to perform an appropriate tuning 
of the couplings in order to cancel the contributions coming from unphysical 
(negative norm) states. 

The inclusion of gauge fields can change the situation |26| , but it is worth- 
while as a probe to see what happens in this limit when negative couplings 
are included, postponing for a future study the effect of gauge fields. 

In this work we study how the existence of opposite couplings influence 
the vacuum of the theory, specifically, whether or not the Ground state (Q) is 
frustrated (the energy cannot be minimized simultaneously for all couplings) 
or even disordered (non zero vacuum entropy). 



2 The Model 

Our starting point is the non-linear a model, with action: 

5 ff = -^#A 



- T+/1 ■ 



(1) 



r,/u 



Where $ is a 4-component vector with fixed modulus <fr r • <fr r = 1. 
The naive way to introduce AF in the non-linear a model is to consider 
a negative coupling. In this case the state with minimal energy for large (3 is 
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a staggered vacuum. On a hypercubic lattice, if we denote the coordinates 
of site r as (r x ,r y ,r z ,r t ), making the transformation 

$ r -> (-l) r -+ r v+ r '+ r ^ T , (2) 

the system with negative (3 is mapped onto the positive f3 one, both regions 
being exactly equivalent. 

Therefore to consider true AF we must take into account either different 
geometries or more couplings, in order to break the symmetry under the 
transformation (0). In four dimensions the simplest option is to add more 
couplings, we have chosen to add a coupling between points at a distance of 
y/2 lattice units. 

Following this we will consider a system of spins {<fr r } taking values in 
the hyper-sphere S 3 C R 4 and placed in the nodes of a cubic lattice. The 
interaction is defined by the action 

S=-P 1 '%2 $ r $ r+/i -fc $ r$r+A+£> , (3) 

r,fi r,n<v 

The transformation (Q) maps the semi-plane (3\ > onto the (3\ < 0, and 
therefore only the region with (3\ > will be considered. On the line /?i = 
the system decouples in two F 4 independent sublattices. 

When 02 — the model is known to present a continuous transition be- 
tween a disordered phase, where 0(4) symmetry is exact, to an ordered phase 
where the 0(4) symmetry is spontaneously broken to 0(3). This transition 
is second order, being the critical exponents those of MFT: a — 0, v — 0.5, 
(3 = 0.5, rj = and 7 = 1 up to logarithmic corrections. The critical coupling 
for this case can be studied analytically by an expansion in powers of the 
coordination number (q = 2d), being f3 c = 0.6055 + 0(q~ 2d ) J27 . 

From a Mean Field analysis, we observe that for /? 2 > the behavior 
of the system will not change qualitatively from the /3 2 = case but with 
higher coordination number. In fact, taking into account that the energy (for 
non-frustrated systems) is approximately proportional to the coordination 
number, there will be a transition phase line whose approximate equation is 

& c + q% = F , (4) 

where q is the quotient between the number of second and first neighbors. 
This line can be thought as a prolongation of the critical point at /3 2 = 
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so the transitions on this line are expected to be second order with MFT 
exponents. This is also the behavior of the two couplings Ising model in this 
region |20 . 

When p 2 < 0, the presence of two couplings with opposite sign makes 
frustration to appear, and very different vacua are possible. 



3 Observables and order parameters 

We define the energy associated to each coupling: 

* = (5 ) 



dp 
0(3, 



1 i-.// 



15^^= (6) 



In terms of these energies, the action reads 

S = -f3 1 E 1 - p 2 E 2 . (7) 
It is useful to define the energies per bound as 

11 

ei = — Ei, e 2 = E 2 , 8) 

where V = L 4 is the lattice volume. With this normalization e\ , e 2 belong 
to the interval [—1,1]. 

We have computed the configurations which minimize the energy for sev- 
eral asymptotic values of the parameters. We have only considered configura- 
tions with periodicity two. More complex structures have not been observed 
in our simulations. 

Considering only the fli > case, we have found the following regions: 

1. Paramagnetic (PM) phase or disordered phase, for small absolute val- 
ues of Pi, (3 2 . 

2. Ferromagnetic (FM) phase. It appears when Pi + Qp 2 is large and 
positive. 

When the fluctuations go to zero, the vacuum takes the form <3? r = v, 
where v is an arbitrary element of the hyper-sphere. 
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Concerning the definition of the order parameter let us remark that 
because of tunneling phenomena in finite lattice we are forced to use 
pseudo-order parameters for practical purposes. Such quantities behave 
as true order parameters only in the thermodynamical limit. In the FM 
phase, we define the standard (normalized) magnetization as 

Mf = ^E # - (9) 

r 

and we use as pseudo-order parameter the square root of the norm of 
the magnetization vector 

m f = (\/m|). (io) 

This quantity has the drawback of being non-zero in the symmetric 
phase but it presents corrections to the bulk behavior order l/\/V. 

3. Hyper-Plane Antiferromagnetic phase (HPAF). It corresponds to large 
Pi, with j3 2 in a narrow interval ([ — /3i/2, — /?i/6] in the Mean Field 
approximation). In this region the vacuum correspond to spins aligned 
in three directions but anti-aligned in the fourth (yu). 

In absence of fluctuations the associated vacuum would be <fr r = (— l) rM v, 
where fi can be any direction, and v any vector on S 4 . 

We define an ad hoc order parameter for this phase as 

MhpAF,, = ^(-l) r ^ r • (11) 
r 

Mhpaf,^ will be different from zero only in the HPAF phase, where the 
system becomes antiferromagnetic on the /i direction. From the four 
order parameters (one for every possible value of fi) only one of them 
will be different from zero in the HPAF phase. So, we define as the 
pseudo order parameter: 



M HP af= .E M hpaf,„- (12) 
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4. Plane Anti-Ferromagnetic (PAF) phase for f3 2 large and negative. In 
this region the ground state is a configuration with spins aligned in two 
directions and anti-aligned in the remaining two. It is characterized 
with by one of the six combinations of two different directions (/i, u), 
and an arbitrary spin v: <fr r = (— l) rM+TV v. For the PAF region we first 
define 

M PAF ^ = i^(-iy> + -$ r , (13) 

r 

and the quantity we measure is 

M PAF= /E M PAF,(^) (14) 



In order to avoid undesirable (frustrating) boundary effects for ordered 
phases, we work with even lattice side L as periodic boundary conditions are 
imposed. 

From this data we can compute the derivatives of any observable with 
respect to the couplings as the connected correlation function with the ener- 
gies 

H = (0£,>-<0><^> (15) 
An efficient method to determine f3 c for a second order transition is to 



measure the Binder cumulant for various lattice size and to locate the 
cross point in the space of (3. 

For O(iV) models U L ((3) takes the form @: 

U L {p) = l + 2/N-^—U (16) 

(irr) 

where m is an order parameter for the transition. 

It can be shown ||, || that U L (0) -> 0(1 /V) and U L (oo) -> 2/N. The 
slope of Ul(P) at (3 C increases with L. 

The value of the Binder cumulant is closely related with the triviality of 
the theory since the renormalized coupling (in the massless thermodynamical 
limit) at zero momentum can be written as: 

g R = lim g n (L) = lim (L/^ L ) d U L ((3 c ) (17) 

L—>oo L—>oo 
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where £l is the correlation length in the size L lattice. 

From this point of view triviality is equivalent to have a vanishing 
in the thermodynamical limit. In this context it is clear that we can use 
the value of g#, to classify the universality class. Out of the upper criti- 
cal dimension, L/£z, is a constant at (3 C since £ ~ L, and we could use the 
Binder cumulant for the same purpose |3(J . At the upper critical dimension, 
£l presents logarithmic corrections and L/£z, is no longer a constant at f3 c . 
For the FM 0(4) model in d = 4 (upper critical dimension) we have per- 
turbatively ~ (InL) -1 / 4 [[nj. In order to have a non trivial theory, 

the Binder cumulant should behave as a positive power of InL, but from its 
definition [|^] we see that Ul(/3) < 1. This is just another way of stating the 
perturbative triviality of the FM 0(4) model. 



3.1 Symmetries on the F 4 lattice 

In the Pi = case the system decouples in two independent lattices, each one 
constituted by the first neighbors of the other. So we consider two lattices 
with F 4 geometry. There are several reasons to choose the point Pi = 
for a careful study of the PM-PAF transition. The region with f3\ > 1.5 
evolve painfully with our local algorithms; For small (3\ we expect very large 
correlation in MC time because the interaction between both sublattices is 
very small, and the response of one lattice to changes in the other is very 
slow. We also remark that the presence of two almost decoupled lattices is 
rather unphysical. 

We also have the experience from a previous work for the Ising model 



20| that the correlation length at its first order transition is smaller in the 
F 4 lattice, that means, we can find asymptotic critical behavior in smaller 
lattices. 

However we should point out that the results in the F4 lattice cannot be 
easily extrapolated to a neighborhood of the Pi axis. Certainly, the geometry 
of the model is very modified when (3\ 7^ 0, and perhaps continuity arguments 
present problems. Nevertheless, we have run also the case Pi ~ 0, and as 
occurs in the Ising model we have not found qualitative differences. 

In the following when we refer to the size of the lattice L on the F4 lattice 
we mean a lattice with L 4 /2 sites. 

We have to find the configurations that maximize E2 in order to define 
appropriate order parameters for the phase transition. 

The system has a very complex structure. As starting point we have 
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studied numerically the vacuum with (3 2 *C 0. For this values we have found 
in the simulation: 

1. The vacuum has periodicity two. To check this, we have defined: 

V * = 7^X>- ( 18 ) 

where % = 0, . . . , 7 stands for the i th vertex of each 2 4 hypercube be- 
longing to the F 4 lattice, and with I we denote the 2 4 hypercubes 
themselves. 

From these vectors we can define the 8 magnetizations associated to 
the elementary cell, 

Vi = (^/V?> , (19) 

We have checked that all V* tends to 1 for the ordered phase in the 
thermodynamical limit, so we conclude that the ordered vacua have 
periodicity two. 

Let us remark for the sake of completeness that all order parameters 
we have defined can be written as an appropriate linear combination 
of the Vi. 

2. In the elementary cell, <fr r+ ^ + p = <fr r V/i, v with /i < v. So, in this 
section we will restrict the study of the vacuum structure to the four 
sites (i — 0, 1, 2, 3) belonging to the cube in the hyper-plane r t = 0. 

3. We have measured the energy per bound associated to the second neigh- 
bors coupling. We check that in the thermodynamical limit e 2 = —1/3. 

4. If we choose the symmetry breaking direction by keeping fix one vector, 
(eg. <&o) we find: 

3 

£((*o-*i)*o-*i) = 0, (20) 

i=i 

The vacuum structure is not completely fixed by these three conditions 
since different symmetry breaking patterns are possible. For instance, a 
configuration $ = (1,0,0,0), $ x = (-1/3, ^ Vl ), $ = (-1/3, ^v 2 ), 
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$0 = ( — 1 /3, ^pv 3 ), with Vi a 3-component unitary vector with the con- 
straint X^j v i v j = 0) breaks 0(4), but an 0(2) symmetry remains (for the 
different Vi). 

To determine which is the vacuum in presence of fluctuations, we consider 
four independent fields in a 2 4 cell with periodic boundary conditions. Let us 
first consider an 0(2) group. We can study the four vectors as a mechanical 
system of mass-less links of length unity, rotating in a plane around the same 
point, whose extremes are attached with a spring of natural length zero. The 
energy for the system is: 

3 

E = - cos ( 9 i - 9 i) ■ ( 21 ) 

i,j=0,i>j 

We consider the fluctuation matrix, H = dE 2 / dOidOj in order to find the 
normal modes. The matrix elements of H take the form: 

H ij = s a cos ( 9i - 9k ) - cos (^ - w - s a) ' ( 22 ) 

In the FM case the minimum correspond to 9{ = 0, for all %. There 
is a single zero mode, and a three times degenerated non-zero mode with 
eigenvalue A = — 4. 

For the AF (maximum energy) case, the maximum energy is found, up 
to permutations, at 9 = <p, 9\ = <p + 7r, 6 2 = <f) + a an d 3 = <p + ir + a, Va. 
In addition to the <fi freedom that corresponds to the global 0(2) symmetry, 
there is a degeneration of the vacuum in the a angle and this zero mode is 
double Va. 

The other two eigenvalues are Ai^ = 2(1 ± cos a), so, an additional zero 
mode appears when a — 0, obtaining in this three fold degenerated 

zero mode corresponding to: 9 = 9\ = 9 2 + n = 9 3 + ir. 

The 0(4) case is qualitatively similar. We have 12 degrees of freedom. Of 
all configurations that minimize the energy, that with a largest degeneration 
(9-times) consist of 2 spins aligned and 2 anti-aligned that correspond to a 
PAF vacuum. We consider this degeneration as the main difference with the 
FM sector, and could be relevant to obtain different critical exponents. 

In presence of fluctuations the configurations with largest degeneration 
are favored by phase space considerations, so we expect that the real vacuum 
is a PAF one. This statement will be checked below with Monte Carlo data 
in the critical region. 
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4 Finite Size Scaling analysis 



Our measures of critical exponents are based on the FSS ansatz [33|. Let 
be (0(L,/3)) the mean value of an observable measured on a size L lattice 
at a coupling (3. If 0(oo,/3) ~ |/3 — /3 c r°, from the FSS ansatz one readily 
obtains |33| 

(0(L,{3)) = L x °/' / Fo{L/aoo,(3)) + ... , (23) 

where Fo is a smooth function and the dots stand for corrections to scaling 
terms. 

To obtain v we apply equation (E3T) to the operator d log M P af / d/3 whose 



related x exponent is 1. As this operator is almost constant in the critical 
region, we just measure at the extrapolated critical point or any definition 
of the apparent critical point in a finite lattice, the difference being small 
corrections-to-scaling terms. 

For the magnetic critical exponents the situation is more involved as the 
slope of the magnetization or the unconnected susceptibility is very large at 
the critical point. 

We proceed as follows (see refs. |T^| for other applications of this method). 



Let be G any operator with scaling law xq = 1 (for instance the Binder 
parameter or a correlation length defined in a finite lattice divided by L). 
Applying eq. ( |2"3"D to an arbitrary operator, O, and to 6 we can write 

(0(L, (3)} = L*°/»f ,e({®(L, P))) + ■ ■ ■ ■ (24) 

Measuring the operator O in a pair of lattices of sizes L and sL at a coupling 
where the mean value of G is the same, one readily obtains 



(0(sL,P)) 



(0(L,(3)) 



= s Xo/u + ... . (25) 

e(L,/3)=6(sL,/3) 



The use of the spectral density method (SDM) |34| avoids an exact a priori 
knowledge of the coupling where the mean values of G cross. We remark 
that usually the main source of statistical error in the measures of magnetic 
exponents is the error in the determination of the coupling where to measure. 
However, using eq. ( plf) we can take into account the correlation between the 
measures of the observable and the measure of the coupling where the cross 
occurs. This allows to reduce the statistical error in an order of magnitude. 
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4.1 FSS at the upper critical dimension: logarithmic 
corrections 



It is well known that d — 4 is the upper critical dimension of the FM 0(4) 
model. As we have already pointed out logarithmic corrections to the Mean 
Field predictions are expected. In particular, FSS in its standard formulation 
breaks at d = 4 because the essential assumption, namely £l(/3 c ) ~ is no 
longer true. In fact, in four dimensions |J: 

£(oo,t) ~ |t|~ 1/2 | ln|t||2 . (26) 

The FSS formula for the correlation length was calculated by Brezin |3l] . 
At the critical point one gets: 

£(L,Pc) ~ L(\nL) 1/4: . (27) 

It has been suggested || that the usual FSS statement should be replaced 
by the more general one: 



0(L,p c ) r /£(£,# 



0(oc,P) ^°V£(oo,/3) 



(28) 



When applying the quotient method described above to systems in four 
dimensions one has to take into account the logarithmic corrections, so that 
the modified formula reads: 



(0(sL,f3)) __/.. A Ins 



(0(L,P)) 



sXO/u ( 1 + ^i) 1/4 ■ ( 29 ) 



e(L,/3)=6(sL,/3) 



This point is particularly important when measuring the magnetic critical 
exponents because as we have already mentioned, the slope of the magneti- 
zation and susceptibility are very large, and one has to be very careful when 
locating the coupling where to measure. 



5 Numerical Method 

We have simulated the model in a L 4 lattice with periodic boundary con- 
ditions. The biggest lattice size has been L = 24. For the update we have 
employed a combination of Heat-Bath and Over- relaxation algorithms (10 
Over- relax sweeps followed by a Heat-Bath sweep). 
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The dynamic exponent z we obtain is near 1. Cluster type algorithms are 
not expected to improve this z value. In systems with competing interactions 
the cluster size average is a great fraction of the whole system, loosing the 
efficacy they show for ferromagnetic spin systems. 

We have used for the simulations ALPHA processor based machines. The 
total computer time employed has been the equivalent of two years of ALPHA 
AXP3000. We measure every 10 sweeps and store the individual measures to 
extrapolate in a neighborhood of the simulation coupling by using the SDM. 

In the F 4 case, we have run about 2 x 10 5 r for each lattice size, being 
r the largest integrated autocorrelation time measured, that corresponds to 
Mpaf, and ranges from 2.3 measures for L = 6 to 8.9 for L = 24. We have 
discarded more than 10 2 r for thermalization. The errors have been estimated 
with the jack-knife method. 

6 Results and Measures 
6.1 Phase Diagram 

We have studied the phase diagram of the model using a L = 8 lattice. We 
have done a sweep along the parameter space of several thousands of itera- 
tions, finding the transition lines shown in Figure 1. The symbols represent 
the coupling values where a peak in the order parameter derivative appears. 

The line FM-PM has a clear second order behavior. It contains the critical 
point for the 0(4) model with first neighbor couplings {f3\ ~ 0.6, = 0) 
with classical exponents [y = 0.5, 77 = 0) . In the /3i — axis, we have 
computed the critical coupling fa 0.18) and the critical exponents as a 
test for the method in the F 4 lattice. We have also considered the influence 
of the logarithmic corrections when computing the exponents. 

The lines FM-HPAF, HPAF-PAF and PM-HPAF show clear metastabil- 
ity, indicating a first order transition. 

The regions between the lower dotted line and the PAF transition line, 
and between the upper dotted line and the FM transition line, are disor- 
dered up to our numerical precision. We could expect always a PM region 
separating the different ordered phases, however, from a MC simulation it is 
not possible to give a conclusive answer since the width of the hypothetical 
PM region decreases when increasing fli, and for a fixed lattice size there is 
a practical limit in the precision of the measures of critical values. 
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On the line PM-PAF we have found no signs of first order. We have done 
hysteresis cycles in several points and no metastability has been observed. In 
Figure 2 we plot the energy distribution at the coupling where a peak in the 
specific heat appears. There is no evidence of two-state signal up to L = 24. 

The likely second order behavior of the PM-PAF transition line contrast 
with the first order one found in the Ising model with two couplings in the 



analogous region pOfl . This is not surprising because we are dealing now 
with a global continuous symmetry. The spontaneous symmetry breaking 
of such symmetries manifest in the appearance of soft modes or low energy 
excitations (long wavelength), the Goldstone bosons in QFT terminology 



35|. The role of these soft modes is quite important and is actually under a 



a vigorous discussion in the two dimensional case [[36], m general, these 



low energy modes will perturb the mechanism of long distance ordering, 
softening in this way the phase transitions. 

Regarding the differences with the FM case, the most remarkable feature 
is the different vacuum structures appearing, specially the very large degen- 
eration in the PAF transition, in contrast with the single degeneration of the 



14 



FM 0(4) mode. 

As the simpler point for study the properties of the transition, namely 
the critical exponents, is the F4 limit, most of the MC work has been done 
for this case. 



6.2 Results on the F 4 lattice 



6.2.1 Results on the FM region 

Firstly, we have checked our method on the FM region of the F4 lattice. In 
Figure 3 the crossing points of the Binder cumulant for various lattice sizes 
are displayed. The prediction for the critical coupling (3 C ~ 0.1831(1) agrees 
with an earlier study by Bhanot [fjSfl . 

Concerning the measures of critical exponents, we have applied the quo- 
tient method, described in section IV. In table [l] we quote the results when 
logarithmic corrections are included (formula (|29|)), and also for sake of com- 
parison, when they are neglected (formula (|25|)). We see how in fact the 
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Figure 2: Energy distribution for L=16,20 and 24 at the peak of the specific 
heat on the F4 lattice. 
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agreement of the critical exponents with the MFT predictions is better when 
the logarithmic corrections are taken into account. 

From now on we will focus on the transition between the PM phase and 
the PAF phase on the F4 lattice. 



6.2.2 Vacuum symmetries on the PAF region 

We will check using MC data that the ordered vacuum in the critical region 
is of type PAF. 
Let us define 

Aij = Vi ■ Vj . (30) 

The leading ordering corresponds to the eigenvector associated to the maxi- 
mum eigenvalue of the matrix A, that should scale as L~ 2l3 ^ u at the critical 
point. The scaling law of the biggest eigenvalue agrees with the (3/v value 
reported in Table ||, and the associated eigenvector is, within errors, (1,1,-1,- 
1). 
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Figure 3: Crossing points of the Binder Cumulant for various lattice sizes on 
the FM-PM phase transition 
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L values 8/16 12/16 10/12 


(without logarithmic corrections) 


ajv 


0.08(5) 


0.02(2) 


0.13(12) 


(5/v 


0.92(3) 


0.94(3) 


0.87(4) 


l/v 


2.16(2) 


2.12(2) 


2.24(4) 


(with logarithmic correc 


tions) 


ajv 


0.0 


0.0 


0.03(8) 


(3/u 


1.04(3) 


1.06(3) 


1.04(2) 


l/v 


1.94(3) 


1.90(4) 


1.93(3) 



Table 1: Critical exponents for the FM-PM phase transition in the F 4 lattice. 



We also have found that the other eigenvalues scale as L~ 4 . This is the 
expected behavior if just the 0(4) symmetry is broken, and it remains an 
0(3) symmetry in the subspace orthogonal to the 0(4) breaking direction. 



6.2.3 Critical Coupling 

To obtain a precise determination of the critical point, /3 C , we have used the 
data for the Binder parameter (p!6|). 

In Figure 4 we plot the crossing points of the Binder cumulants for the 
simulated lattices sizes. Extrapolations have been done using SDM from 
simulations at (3 2 = -0.7090 for L = 6,8, 10, 12 and 16; (3 2 = -0.7078 for 
L = 20, and (3 2 = -0.7070 for L = 24. 

The shift of the crossing point of the curves can be explained through 
the finite-size confluent corrections. The dependence in the deviation of the 



crossing point for L and sL size lattices was estimated by Binder |28 



(3 C (L, sL) - (3 C ~ ^-f-^L—i/" , (31) 

where u is the universal exponent for the corrections-to-scaling. 
The infinite volume critical point the value 

/3 C = -0.7065(5) [+2] [-2] , (32) 

where the errors in brackets correspond to the variations in the extrapolation 
when we use the values u = 0.5 and u = 2 respectively. In Figure 5 we plot 



eq. (|3lD for u = 1. 
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0.50 r-r 
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-0.710 -0.708 -0.706 

Figure 4: Crossing points of the Binder cumulant for various lattice sizes. 

Using the previous value of j3 c we can compute the Binder cumulant at 
this point. In table § we quote the obtained values. The result points to that 
the Binder cumulant stays constant in the critical region. This result would 
be compatible with a non zero value of the renormalized coupling when L 
increases. 

Concerning the possibility of having logarithmic corrections in the deter- 
mination of the critical coupling, from the numerical point of view, it is not 
possible to discern between the u effect, and a logarithmic correction. 



6.2.4 Thermal Critical exponents: a, v 

The critical exponent associated to correlation length can be obtained from 
the scaling of: 

« = -^-. (33) 

where M is an order parameter for the transition, Mpaf for our purposes. 
In the critical region k ~ L l / U . As k is a flat function of /3, is not crucial 
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the point where we actually measure. The results displayed in table |3| have 
been obtained measuring at the crossing point of the Binder parameters for 
lattice sizes L and 2L using (p5|) . 

For measuring a/v we study the scaling of the specific heat 

We expect that C scales as A + BL a ^ u , where A is usually non- negligible. In 
Figure 6 we plot the specific heat measuring at (p2|), as well as at the peak 
of the specific heat, as a function of L. We observe a linear behavior for 
intermediate lattices. For the largest lattice the slope decreases. The weak 
first order behavior (a/v = 1 for small lattices that becomes d for large 
enough sizes) seem hardly compatible with our data. If we neglect the A term 
(what is asymptotically correct), and compute the exponent using eq. ( |25|) 
we obtain a/v ~ 0.3 for intermediate lattices that reduces to a/v = 0.15(2) 
for the (20,24) pair. 

However, to give a conclusive answer for the value of a statistics on larger 
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Lattice sizes 


U L (p c (u = 0.5)) 


U L (p c (u=l)) 


U L (p c (u = 2)) 


6 


0.4435(15) 


0.4437(12) 


0.4438(12) 


8 


0.4406(15) 


0.4409(12) 


0.4413(12) 


10 


0.4407(14) 


0.4411(16) 


0.4414(14) 


12 


0.436(4) 


0.437(3) 


0.438(3) 


16 


0.435(3) 


0.436(3) 


0.437(3) 


20 


0.429(5) 


0.431(5) 


0.433(5) 


24 


0.428(6) 


0.430(7) 


0.433(7) 



Table 2: Binder cumulant for various lattices sizes at the extrapolated 
critical point for u = 0.5, 1, 2. 
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Figure 6: Specific heat at the peak (triangle symbols) and at P — —0.7068 
(cross symbols) as a function of the lattice size. 



lattices are mandatory. 
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0.98 0.99 1.00 1.01 1.02 




0.98 0.99 1.00 1.01 1.02 



U 2i/Ul 

Figure 7: Quotients to obtain (3/v and j/u. 



6.2.5 Magnetic Critical exponents: 7, (3 

The exponents 7 and (3 can be obtained respectively from the scaling of 
susceptibility and magnetization: 

X = V(M 2 ) ~ V lu (35) 

M ~ L"^" (36) 

Where M is an order parameter for the phase transition. In Figure 7 
upper part, we plot the quotient between M P af for lattices L and 2L as a 
function of the quotient between the Binder cumulants for both lattice sizes. 

For large L in the critical region we should obtain a single curve, the 
deviations corresponding to corrections to scaling. In the lower part of Figure 
7 we plot the same function for susceptibility. 

The values for 7 and (3 are summarized in Table |3|. 
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Lattice sizes 7/^ 


P/v 


(without logarithmic corrections) 


6/12 


2.417(3) 


0.791(4) 


0.474(10) 


8/16 


2.403(3) 


0.792(6) 


0.483(8) 


10/20 


2.410(2) 


0.790(4) 


0.471(6) 


12/24 


2.403(5) 


0.797(5) 


0.483(7) 


20/24 


2.398(5) 


0.802(4) 


0.487(6) 


(with logarithmic corrections) 


6/12 


2.301(3) 


0.849(4) 


0.484(9) 


8/16 


2.300(3) 


0.850(5) 


0.489(7) 


10/20 


2.315(2) 


0.843(3) 


0.488(5) 


12/24 


2.314(2) 


0.842(5) 


0.487(5) 


20/24 


2.317(5) 


0.839(4) 


0.498(5) 



Table 3: Critical Exponents for the PM-PAF phase transition in the F 4 
lattice 



6.3 Logarithmic corrections 

We now address the question of the possibility of logarithmic corrections in 
the AF 0(4) model. For the thermal critical exponents, the situation seems 
clear, they are compatible with the classical exponent v = 0.5. For the 
magnetic exponents, the situation is more involved. In principle, one can 
think that they disagree from MFT due to logarithmic corrections. We have 
no perturbative predictions about the form in which these corrections would 
affect £l for the AF case. However, one expects that such corrections slightly 
modify the critical exponents, as occurs in the FM case. It could be possible 
that logarithmic corrections modify largely the previous critical exponents 
and drift them to the FM ones. To sort this out, we have considered the 
possibility of a behavior FM like, so that £l ~ L(lnL) 1 / 4 . In the lower part 
of Table we quote the values of the critical exponents for the PAF phase 
transition when logarithmic corrections are included (formula (^)). We see 
how in effect the magnetic critical exponents are too far from the classical 
ones for being the result of a logarithmic correction to the MFT predictions. 
It is interesting to compare this situation with that in the RP 2 model in 



d = 4 |14| where small deviations from MFT exponents can be explained as 



logarithmic corrections. 
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7 Conclusions and outlook 



We have studied the phase diagram of the four dimensional 0(4) model with 
first and second neighbors couplings. For < we find a region non-trivially 
related with the FM one, in which the system is AF ordered in some plane. 
The phase transition between the disordered region and this PAF region 
seems to be second order. 

We also compute the critical exponents on this line at /?i = (F4 lattice) 
by means of FSS techniques. We found that up to L = 24 the exponents are 
in disagreement with the Mean Field predictions. Specifically, from our 7/1/ 
estimation (or /3/v using hyper-scaling relation) the exponent 77 associated 
with the anomalous dimension of the field is 77 ~ —0.4. This fact itself 
would imply the non-triviality of the theory because Green functions would 
not factorize anymore. One cannot discard that the observed behavior were 
transitory. However, the stability of our measure of 7/V for lattice sizes 
ranging from L = 6 to L = 24, which are more than a hundred of standard 
deviations apart the MF value, makes very unlikely this hypothesis. Actually, 
it would be possible to obtain triviality also with a logarithmic exponent in 
equation (|29| ) different from 1/4. We can fix the critical exponents to its 
MF value and compute this parameter from the numerical data. The results 
obtained shown a non asymptotic behavior, with values ranging from 0.8 
to 1.2 for the lattices used. A logarithmic fit is not satisfactory because of 
the non-asymptoticity and the large value of the logarithmic exponent, but 
larger lattices sizes are needed in order to get a more conclusive answer. 

The behavior of the specific heat does not show any first order signature, 
but we have been not able to obtain a reliable estimation of the a exponent. 

We have also measured the Binder cumulant at the critical point, finding 
that it stays almost constant when increasing the lattice size. If this is not 
a transient effect, and logarithmic corrections are finally ruled out, it would 
correspond a nonzero value of the renormalized constant in the thermody- 
namical limit. 

8 Acknowledgments 

Authors are grateful to J.L. Alonso, H.G. Ballesteros, V. Martin-Mayor, J.J. 
Ruiz-Lorenzo, CL. Ullod and D. Iniguez for helpful discussions and com- 
ments. We thank specially J.M. Carmona for useful discussions concerning 



23 



the logarithmic corrections. One of us (I.C.) wishes to thank R.D. Kenway 
his kind hospitality at the Department of Physics and Astronomy, Univ. of 
Edinburgh, as well as to Edinburgh Parallel Computing Center (EPCC) for 
computing facilities where part of these computations have been done under 
the financial support of TRACS-EC program. This work is partially sup- 
ported by CICyT AEN94-0218, and AEN95-1284-E. I.C. holds a Fellow from 
Ministerio de Education y Ciencia. 

References 

[1] D. Brydges, J. Frohlich and T. Spencer Commun. Math Phys. 83, p 
123 (1982) 



[2] 
[3] 
[4 



[5] 
[6] 
[7 
[8 
[9 
[10 

[11 
[12 



M. Aizenman Phys. Rev. Lett. 97, p 1 (1981) 
J. Frolich, Nuc. Phys. B200 [FS4], p 281 (1982) 

E. Brezin, J.C. Le Guillou and J. Zinn- Justin. "Field Theoretical ap- 
proach to critical phenomena" in Phase Transitions and Critical Phe- 
nomena ed. C. Domb and M.S. Green (Academic Press, London) (1976). 

F. J. Wegner and E.K. Riedel Phys. Rev.B7, p 248 (1973) 

M. Luscher and P. Weisz Nuc. Phys. B299 [FS20], p 25 (1987) 
D. E. Callaway, Nuc. Phys. B233, p 189 (1984) 

RTN Collaboration: J.L. Alonso et al, Nuc. Phys. B405, p 574 (1993) 

R. Kenna and C.B. Lang Nuc. Phys. B393, p 461 (1993) 

R. Kenna and C.B. Lang Nuc. Phys. B30 p 697 (Proc. Suppl.) (1993) 

I. Montvay Nuc. Phys. B26 p 57 (Proc. Suppl.) (1991) 

H. G. Ballesteros, L.A. Fernandez, V. Mart in- Mayor and A. Munoz- 
Sudupe, Phys. Lett. B378, p 207 (1996); \hep-lat/960503^ to appear 
in Nuc. Phys. B 



[13] E. Seiler and K. Yildirim hep-lat ^-preprint 96090~3(\ 



24 



[14] H.G. Ballesteros, J.M. Carmona, L.A. Fernandez, V. Martin-Mayor 
A. Muiioz Sudupe and A. Tarancon. preprint DFTUZ 96/20 

[15] M.L. Plumer and A. Caille, J. Appl. Phys. 70, p 5961 (1991). 

[16] H. Kawamura, J. Phys. Soc. Jpn. 61, p 1299 (1992). 

[17] H. Mori and M. Hamada, Physica B 194-196, p 1445 (1994). 

[18] J.L. Moran-Lopez, F. Aguilera-Granja and J.M. Sanchez, J. Phys.: 
Cond. Matt. 6, p 9759 (1994). 

[19] G. Parisi and J. J. Ruiz-Lorenzo, J. Phys. A (Math and Gen.) 28, p 
L395 (1995). cond-mat preprint 9503016. To appear in J. Phys. A. 

[20] J.L. Alonso, J.M. Carmona, J. Clemente, L.A. Fernandez, D. Ihiguez, 
A. Tarancon and C.L. Ullod, Phys. Lett. B376, p 148 (1996) 

[21] W. Bock, A.K. De, K. Jansen, J. Jersak, T. Neuhaus and J. Smit Nuc. 
Phys. B344, p 207 (1990) 

[22] W. Bock, A.K. De, C. Frick, J. Jersak and T. Trapenberg Nuc. Phys. 
B378, p 652 (1992) 

[23] T. Ebihara and K. Kondo Prog. Theor. Phys. vol 87, 4 (1992) 

[24] J.L. Alonso, Ph. Boucaud, F. Lesmes and E. Rivas. Phys. Lett. B329, 
p 75 (1994) 

[25] G. Gallavotti and V. Rivasseau, Phys. Lett. B122, p 268 (1983) 
[26] Jochen Fingberg and Janos Polonyi, hep-lat preprint 9602003 
[27] P.R.Gerber and M.E. Fisher, Phys. Rev. BIO, p 4697 (1974) 
[28] K. Binder, Z. Phys. B43, p 119 (1981) 

[29] E. Brezin and J. Zinn- Justin Nuc. Phys. B257 [FS14], p 867 (1985) 
[30] Giorgio Parisi, Statistical Field Theory, Ed. Addison Wesley (1992). 
[31] E. Brezin J. Phys. 43, p 15 (1982) 



25 



[32] C. Itzykson and J.M.Drouffe, Statistical Field Theory, Vol. 1, Cam- 
bridge University Press (1989) 

[33] Finite size scaling, Current Physics - Sources and comments, vol. 2, Ed. 
J.L. Cardy, Elsever Science Publishers B. V. (1988) 

[34] A.M. Ferrenberg and R.H. Swendsen, Phys. Rev. Lett. B61, p 2635 
(1988) 

[35] Michel Le Bellac Des phenomenes critiques aux champs de jauge Inter- 
Editions /Editions du CNRS, Paris (1988) 

[36] B. Alles, A. Buonanno and G. Cella hep-lat ^-preprint 9609030] (To 
appear in the Lattice '96 Proc. Suppl. Nuc. Phys. B) 

[37] E. Seiler and A. Patrasciou hep-lat ^-preprint 9608138] 

[38] G. Bhanot Nuc. Phys. B17, p 653 (Proc. Suppl.) (1990) 

[39] L.A. Fernandez, M.P. Lombardo, J.J Ruiz-Lorenzo and A. Tarancon, 
Phys. Lett. B277, p 485(1992) 



26 



